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We discuss the phase structure (in the 1/A^-expansion) of the Nambu-Jona-Lasinio model 
in curved spacetime with non-trivial topology x S^. The evaluation of the effective po- 
tential of the composite field ipip is presented in the linear curvature approximation (topology 
is treated exactly) and in the leading order of the 1/iV-expansion. The combined influence 
of topology and curvature to the phase transitions is investigated. It is shown, in particular, 
that at zero curvature and for small radius of the torus there is a second order phase tran- 
sition from the chiral symmetric to the chiral non-symmetric phase. When the curvature 
grows and (or) the radius of S^ decreases, then the phase transition is in general of first 
order. The dynamical fermionic mass is also calculated in a number of different situations. 
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1 Introduction 



The study of composite fermionic fields — and effects related with them — in the very early 
universe has attracted much attention in the physical community. In order to be able to 
do such kind of research, one must first develop an effective action formalism for composite 
fields in a general curved spacetime (for an introduction to this subject, see [|l|). However, 
even in flat space only very few models are known which can be treated analytically, when 
studying the composite boundstates. The Nambu-Jona-Lasinio (NJL) model (for a recent 
discussion, see |^) belongs to such interesting class. This model is usually discussed in frames 
of the 1/A^-expansion (see, for example 0]), which is a very useful scheme to study the non- 
perturbative effects. The NJL model may be considered as an effective field theory for QCD 
in some region, and it is also connected with the theory of superconductivity (related models 
with quite interesting properties were studied some time ago p). One can explicitly discuss 
the dynamical chiral symmetry breaking for the NJL model in the 1/A^-expansion, where a 
vacuum condensate < ipip >^ appears and a dynamical fermionic mass is generated. 

Recently [^, 0], a detailed study of the NJL model in curved spacetime has been started. 
The one- loop effective potential in the 1/A^-expansion and in the linear curvature approxima- 
tion has been calculated and the existence of a curvature-induced first-order phase transition 
from a chiral symmetric to a chiral non-symmetric phase has been shown (for a review 
on curvature- induced phase transitions for elementary fields in GUT's see, for example ||1|]). 
However, one might expect that the very early universe had a non-trivial topology and (or) 
that it was very hot. This renders interesting to investigate the phenomenon of dynamical 
chiral symmetry breaking in curved spacetime with non-trivial topology. 

The present work is specifically devoted to the study of the NJL model in a curved 
spacetime of the form Ai^ x S^, where is the one-dimensional sphere and Ai^ a three- 
dimensional, arbitrarily curved manifold of trivial topology. Here, we will only consider the 
standard choice of boundary conditions (periodic and antiperiodic) for the fermion on the 
one-dimensional sphere p[. However, more general choices of the boundary conditions 
are, of course, possible 0, in particular those in which the fermion ip acquires an extra phase 
exp {2TT(f) with arbitrary (f (see the second ref. in 0) each time it goes along S^. The (f = 
case corresponds to periodic boundary conditions and (p = 1/2 corresponds to antiperiodic 
boundary conditions. 

In the next section we calculate the effective potential for the NJL model in the spacetime 
Ai^ X S^, in the 1/A^-expansion and for the linear curvature approximation. In section 3, the 
phase structure of the effective potential is discussed and the dynamically generated fermionic 
mass is calculated. The existence of topology- and curvature-induced phase transitions from 



2 



the chiral symmetric phase to the non-symmetric one is estabhshed for different regions of 
the parameters of the theory. 



2 The effective potential for the NJL model 



Let us start the calculation of the effective potential for the NJL model in the 1/A^-expansion. 
The classical action of the theory in curved spacetime is given by 

A 



S 



2N 



where is the number of fermions and the rest of the notation is standard @, 0. It is 
more useful in actual calculations to work with the following equivalent action, in which the 
auxiliary fields a and tt are introduced: 

N 



S 



2A 



(2) 



Already in flat spacetime is it known that the global chiral symmetry — which is the classical 
symmetry of the theory (Q) or (0) — is spontaneously broken when the coupling constant A 
exceeds some critical value Ac- The contribution of the external gravitational field to this 
effect has been discussed in ref. 0. Our purpose here will be to determine the influence 
of combined effects, namely of external gravity and non-trivial topology, simultaneously, on 
the dynamical chiral symmetry breaking and restoration. Notice that such a study cannot 
be done for the most simple low- dimensional analogue of the NJL model, namely the D = 2 
Gross-Neveu model neither for the — a little more complicated — D = 2 Schwinger |11] 
or D = 2 Thirring |jl2| models. In all those cases one can study only the isolated effects of 



either non-trivial topology (non-zero temperature [|T3|) or external gravity fT^ (see also 0) 
on the chiral symmetry breaking pattern. 

Let us now proceed with the explicit calculations. First of all, we introduce the semiclas- 
sical effective action Sejf as 

Z[0,0] = j Da DnexpiiNSeff) , (3) 

where Z is defined, in the usual way, to be the generating functional 

Z[r], v] = J Dip Dip Da Dn exp (iS + if]ip + iipr]^ . (4) 

This Seff is trivially solved by integrating over the variables ip and ip, what amounts to a 
simple gaussian integration, namely 



r r 1 

Seff = -J d^x^ — (d^ + vr^) + ilogdet [ii^{x)V^ - {a + i-f^n)] . 



(5) 
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Moreover, Seff is the leading term in the large- expansion of the effective action F. 

The effective potential (cr, vr) is defined by (a, vr) = — F (a, vr) /{N ■ Volume), with 
constant configurations for the fields. Thus, to leading order in the 1/A^-expansion, we have 

V (a, tt) = ^ (fT^ + TT^] + i Sp log < X I [i-i^" (x) - (a + 27577)] \ x> (6) 

The second term is directly related to the Green function given by the equation 

[ir{x)V,,-s]S{x,y)=5\x,y), 

where 5 refers to the scalar, coordinate independent Dirac delta functional for the given 
manifold. This fact can be seen by directly applying the operator expression 

log ( — - — \ = — [ dr- 



A ) Jo A-r 
to i). 

To calculate the second term of the Ihs of equation (||), we make use of Schwinger's 
proper-time method. First, we write 



V {a, tt) = ^ [cr' + Sp log < x | {i-f^V^, - s) ' | x >, (7) 

and using the known expression for the propagator of a free Dirac field in a weakly-varying 
gravitational background (which has been obtained in terms of the Riemann normal coor- 
dinate expansion fl^), we will then write the effective potential with accuracy up to linear 
curvature terms. Notice also that a summation over the two inequivalent spin structures 
which are admitted by the spacetime Ai^ x will be performed P]. Of course, one can con- 
sider also the corrections corresponding to periodic and antiperiodic (non-zero temperature) 
boundary conditions independently. We will make a few comments about this point later. 
Thus, the effective potential is found to be 



- — R h^'ka + S) ^ + -Rauk^'k" h^'ka + s) ^ ^ 



---f^r'^R^dauk" — ^ — 
2 ' (A;2 - s2) 

where one should integrate over k^,k^,k'^ and sum over the coordinate fc^, which is given 
by: k^ = {2n + 6p^i) tt/ L. In expression (||), tr only refers to the spinor indices. We have 
set TT = 0, since there is a rotational symmetry in the fields a and tt, so that it is enough to 
discuss the a 7^ case for the effective potential only. 
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Integration over s is immediate. To perform the momentum integration, one first makes 
the Wick rotation = ik"^) and puts then a cut-off to regularize the resulting expressions. 
In our case, we simply restrict 



so that our cut-off is different, when compared with the cut-off for the case of trivial topology 
i- 

From now on, we shall call Vi the contribution to the effective potential which comes 
from the logarithm of the determinant of the operator which appears in equation (^. After 
carrying out the integrations over s and the momenta, we are led to the following expression 
for the contribution to Vi coming from the p = case, namely, purely periodic boundary 
conditions (this corresponds to the contribution to V{(T, 0) obtained by taking only the p = 
term in equation (j^)). Of course, we want to study Vi, which is given by the sum over the 
two values of p, but — as we discuss below — Vi may be written down immediately once the 
p = contribution has been worked out. We obtain 



p=0 



-| n=oo 

I S 

n=— CO 



f 



Svr' 



A^log 1 + 



A2 + ^ 



37r2 
R 



3/2 



L2 



arctan 



A 



L2 + O- 



+ 



3(27r) 
R 



2Tm 



■ arctan 



AL 

2?™ 



cr2 arctan 



A 



arctan 



(—) 

\27mJ 



L2 



+ 0-2 



2(37r) 



A2 + 



0^ + (t2 + A2 



(9) 



To simplify this expression we use standard techniques drawn from complex analysis, 
such as the expression 



^ 27171 \ L 



27rz 



ioo X L 
dp-[fip) + fi-p)] + — 

ioo I ZTTI 



dpMppl. (10) 

: ^ exp {Lp) - 1 ^ ' 



One has just to identify the function / for each term in equation (P) and then perform the 
integrals in (|10D- Another remark is in order here: when computing the term 



2 1 
L3^ 



E 



2^2 



a2A 



3/2 



arctan 



A 



/47rfrif 
L2 



+ a2 



2™ 



arctan 



AL 
27m 



it is better to rewrite it as 



since now the expression within square brackets satisfies the properties which justify the use 
of equation ([T0|). 

As we said before, once Vi~^ has been computed, one may write immediately Vi~^ and 
Vi, which is given (by definition) by Vi = Vi~^ + Vi~^, because 



-| oo 

7 E 



2n + l 



(In 



and, then 



-| oo 

7 E 



2n + 1 



7EW7 



(11) 



(12) 



Now it is apparent that the physics displayed by the model in the case of purely periodic 
boundary conditions and in the one where we consider both spin structures will be essentially 
the same, since 



T/P=0 I yl 



p=l 



2V[, 



2L ■ 



Thus we see that both cases are related by a trivial rescaling of the length and an overall 
factor which multiples Vi. That is quite an appealing result. 

However, the last remark does not apply to the case of purely antiperiodic boundary 
conditions, where we only take p = 1 — which gives the themodynamics of a system in 
three-dimensional space with trivial topology. In fact, from the last expression we get 



Vl,L - Vff 



2 '2 



Henceforth, we shall concentrate below only on the analysis of the case in which both spin 
structures are taken into account, as they appear in equation (|^). 

The first term on the rhs of equation ( p!OD may be computed without difficulty in all cases. 
The second term is, in general, rather more involved. In order to simplify this contribution as 
much as possible in the cases which come from equation (|^), one should pay careful attention 
to the determination of the integrand along the contour of integration. After some work, the 
final result is found to be 



A4 



37r2 
1 



1 - Vl + X2 + - log 



'exp {21^/1 + 



- 1 



4 
3^ 



x^ 



VT- 



exp (2/) 
3 



x^ 



-X arcsinh — 



X 



dr 



{t' - x') 
exp (2/r) 



dr- 



exp (2/r) 
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6 (27r)^ 
2r 



1 — vT+~x^ + x^arcsinh 



3 (27r)' [■ 
r 



Jo exp (2Zr) — 1 Jx 
1 



(13) 



dr- 



2 (37r)^ 



+ 



exp (2lr) - 1 
1 1 



VTT^ exp (2^^/^T^) - 1 exp (21) - 1 



where x = a/ A, I = LA and r = R/ A^. 

The value of the field cr which satisfies the gap equation 

gives a dynamical mass to the fermions. This last equation, when written in terms of the 
natural variables x, I, r and c (c = AA^), reads 

V'(x) x^ 5 X 4: X 1 







A4 



2c 67r2 v/r+^ 37r2 ^1+^ 1 - exp (-2/^1+^) 



X 

+ — 



X arcsinh 



+ 



X 



37r2 

4 



VT+^exp (2lVTTx^) - 1 



27r2^/rqr^ 

— 3a; 



dr 



2\V2 



x^ 



exp(2ZT) - 1 



X 



(14) 



3(27r) 
2rx 

3 {27if [VT+^exp (2/x/rT^) 

,x 1 



4 (37r)' (1 + x2)='/' 



1 Jx 



dr 



Vt2 - x2 exp(2/T) - 1 



2/ 



exp 2/Vl + x'^ 



2 {37rf [{1 + x^f exp (2/^1+^) - 1 I + (exp (2^^!+^) - l)' 

It is not possible to produce an exact, analytical expression for the dynamically generated 
fermion mass. Therefore, we will calculate it below only in a number of limiting cases. 



3 Small-L limit 

One gets convinced immediately that it is very difficult to go any further without doing 
some simplification. In this section we will consider the case LA <C 1, and we will treat the 
opposite case LA ^ 1 in the next one. 



Let us expand V/A'^ in powers of /. Assuming now that + x"^ <C 1, we readily obtain 
IVix) x' 1 . , ,x , 2 / , , I , 1 



^log(l + .V^(- 



A^ 2g Sn^-^y-'^J ' 3.M-" + 
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H ^^a;arcsecWl + ^ H — ^-^ - l] + O (l'^) , (15) 

3(27r)' V 2(37r)' \1 + J W ' ^ ^ 

where (7 = \h? /I. In order to study the phase space of the model, one has to compute also 
the derivative of the effective potential. That yields 



IV' X 2 X 2 ^ / 1 4 X 



2 



H — -X arcsecW 1 H — r 



A*^ g Stt^ 1 + vr^ V x"^ Svr^ 1 + x^ 

r / 1 r X r X . 

H ^-arcsecWlH — ^ tj -. (16) 

3(27r)' V 3^=^ (37r)'(l + x2)' 67r2 1 + x^ ^ ^ 

3.1 Phase structure in the fixed-curvature case 

We will here analyze the phase structure of the model in the case of fixed curvature. In par- 
ticular, the situations of zero curvature and constant non-zero curvature will be considered. 

3.1.1 Zero-curvature case 

As the validity of our results is restricted by the condition r <C 1, it is worthwhile studying 
first the case r = 0. Here it can be proven that it is impossible to have a first order phase 
transition, in fact there is a second order phase transition at g^r = 7i'^/2. For g > g^r there 
is chiral symmetry breaking and the symmetry is restored for g < gcr] in other words, for 
a given value of A, the symmetry is restored when / grows beyond a critical point. In the 
broken phase the order parameter is given by xtr — ~ '^^j terms of the generated 
mass mgen = ^ — -y-. It is straightforward to study the behavior of the value of the effective 
potential at Xbr, for varying g or I, in this limit of small compactification length. One finds 

f K) - -4 (li - i)' . 

This last expression has the appearance of some kind of dynamical dimensional reduction 
(see ref. "^^^^ ^^^^ vacuum becomes less and less energetic as the 

compactification length shrinks to zero (see Fig. 1). In this limit we obtain: 

3.1.2 Influence of the curvature 

Now we shall study the influence that the presence of curvature has on this behavior. It 
is immediate to notice that for negative values of the curvature, the symmetry is always 
broken (the slope of the effective potential is r/(247r) at the origin, see Fig. 2). Moreover, 
for fixed positive values of r (which is kept always small), as I grows there is a first-order 
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phase transition: it has actually lost its continuous character. Now the critical value of the 
parameter g is 

gcr = 



2 V 3 

The approximation is consistent with the small-curvature limit. The difference between the 
values of the order parameter x in the broken and disordered phases at the phase transition 
is given by x = r/8 (retaining again only the first correction coming from the curvature). 
The influence of curvature on chiral symmetry breaking m D = 2 fermionic models is very 
similar (see refs. fl], 0). 

3.2 Phase structure in the case of fixed compactification length 

We can now, on the other hand, study the situation when g is held fixed and the curvature 
takes on different values. 

3.2.1 Case g < n'^ /2 

If g < 7T^/2 there is a continuous phase transition at r = (the symmetry is broken for 
negative values of r and is restored for positive curvature). As r approaches from below, 
the order parameter tends to zero according to the expression 



X 



247rA 



or rrigen = being A = ^ — In deriving this result we assume that | r |^ 1 and 

that I r |< A. 

3.2.2 Case g > 7i^/2 

To finish this section, we consider the case A < or, equivalently, g > 7r^/2. To be consistent 
with the requirement of small curvature near the critical point, we have to assume now that 
I A 1^ 1. In this setting one may expect that, as r grows, there will be a phase transition 
from the ordered to the disordered phase at some positive value of the curvature. Retaining 
again the first terms of the expansions only, we obtain 

(Stt)^ 

and the value of the order parameter at the transition is (for the ordered phase) 

X = — A^ . 

4 
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4 Large- L limit 



We shall here consider again the phase structure in the fixed-curvature case but corresponding 
now to the opposite limit, when L is large. 

4.1 Zero-curvature case 



It is easy to see from equation (|T^) that, dropping exponentially vanishing contributions, 
one may approximate it with 

Vix) If/ , ^ \ 2 / 9 4 1 1" 

— — ;— = t; 2 [vl + — 1] + X yl + x-^ — x arcsinh— , 

where c is, as before, c = AA^. 

Let us consider the gap equation V'{x) = 0. Differentiation of eq. (0) yields 



V'{x) _ X (tx^ 



Vi + a? + x^ arcsinh— 

\ c X ^ 

From here it is trivial to see that the symmetry is broken when 
Otherwise, the symmetry is respected by the vacuum. 



:i7) 



(18) 



4.2 Influence of the curvature 

If one takes into account the presence of a background gravitational field, one can check that 
the terms missing from equations (|17D and (pIS]) are such that now 



V{x) 



X 



2c 47r2 
r 



2 ( vT+a? — 1 ) + x^yl + x^ — x'^arcsinh 



1 



X. 



6(27r2 



1 — vT+lr^ + arcsinh— ^ 

xj 



2(37r) 



x^ 



and 



V\x) 
A4 



X 



a/1 + x^ + x^arcsinh- 

c X 



12 V 



x^ 



18(1 + x 



2^3/2 



(19) 



(20) 
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4.2.1 Case A < ttVA^ 

In this situation one may expect that there will be a continuous phase transition. After 
a short calculation, expanding equation (|20|) around the origin and keeping only the two 
first leading contributions to V'{x), one sees that, very near the phase transition, the order 
parameter is found by solving 

TT^ 5 r , 1 / 2 1 1 ~ 

= 1 r H arcsmh — \- O [x arcsinh — 

c 36 12 X \ X 



Thus, in a first approximation 



X 



sinh(i^(vrVc-l)-^ 



(remember that this is valid for negative r). The symmetry is restored for positive values of 
r. 

4.2.2 Case A > tt^/A^ 

Here, a quick analysis of the previous expressions tells us that in this case there is a first- 
order phase transition, from the disordered to the ordered phase, as the curvature grows 
beyond a positive critical value. Notice that this result is very much like the one that one 
obtains in the small-L limit (see Fig. 3). 



5 Conclusions 

Recently there has been an increased interest in the NJL model, in connection with the dy- 
namical symmetry breaking of the electroweak interactions, using the top quark condensate 
as order parameter [Q, |18|. Here, we have investigated the phase structure of this model in a 
curved spacetime with non-trivial topology — using the 1/A^-expansion and working in the 
linear-curvature approximation — and we have explicitly shown the possibility of curvature- 
and (or) topology-induced phase transitions from a chiral symmetric phase to a chiral non- 
symmetric one. In our approximate analysis, we have discussed the composite field effective 
potential where summation over two inequivalent spin structures has been performed. Of 
course, one can get the results corresponding to purely periodic (or antiperiodic) boundary 
conditions as some particular cases of the above study. 

There are different possibilities to extend our analysis. First, one can consider different 
topologies, for example, Ai^ x T^ or the hyperbolic one Ai^ x TH? /T . Second, it would be 
of interest to study the same problems treating the external gravitational field exactly. In 
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particular, one possibility is to take the De Sitter space S4 for such a calculation. These 
considerations may draw some connection with renormalized quantum gravity in the 
expansion (see [0), that could be certainly important for quantum cosmology near the 
Planck scale. We plan to discuss these questions in the near future. 
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Figure captions 

Figure 1. Plot of the function V/A^ for fixed c = 0.05, r — 0, and different values of I. 

Figure 2. Plot of V/A^ for fixed g — 5.0404 (case g > 7r^/2) and for different values of the 
curvature. It shows the discontinuous character of the phase transition in this case. 

Figure 3. Plot of the effective potential in the large-L limit, corresponding to the case 
A > TT^/A^ (c is held fixed at c = tt^ + 0.15) for different values of r. The discontinuous 
character of the phase transition is clearly seen. 
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